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1. INTRODUCTION 

Now in the word, methods of medical image acquisition methods have advanced rapidly, including 
magnetic resonance imaging (MRI), computed tomography (CT), ultrasound (US), positron emission 
tomography (PET), and single photon emission tomography are some of the imaging techniques. Markov 
random fields in image segmentation is a book that covers the principles of markovian modeling in image 
segmentation and provides a quick summary of current achievements in the subject. Image labeling is a 
common framework in which segmentation is considered, when the task is simplified to labeling pixels. With 
the production more and more of image medical, image automated, processing and analysis techniques have 
become increasingly important [1]. In the field of probability, it is possible to introduce a segmentation class 
where the characteristics of regions can be modeled by MRF. The major drawback of the Markov field model 
is that the approximation algorithms they induce are iterative and very computationally intensive. To reduce 
the amount of computation and obtain non-iterative computational algorithms, we can use hidden Markov 
chains [2]-[4], or with Gaussian mixture models (GMM) [5], [6], or with conditional random fields 
conditional random field (CRF) [7], [8]. 

Although optical, robotics has made extensive use of image processing technology most of these 
methods cannot be used to process medical images. Image segmentation is the process of dividing an image 
into multiple segments to simplify the representation of the image or extract meaningful objectsImage 
segmentation is one of the most important roles in computer vision, and it has several applications in domains 
such as pattern recognition, remote sensing, and machine learning, medical diagnostics and computer vision, 
remote sensing, and medical imaging. A fuzzy treatment of masked image segmentation with MRF model 
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has been proposed [9]-[12]. Thus, this work is novel in three ways. i) we offer an estimation and prediction of 
interface levels. Approach based on MRF, ii) a method for simultaneously estimating GMM parameters using 
the EM Algorithm to account for MRF spatial constraints; and iii) validation of proposed methods on medical 
and synthetic images [13]. 

The MRF model MRF [14] is a probabilistic graphical model that provides a statistical framework 
for modeling spatial contextual constraints as prior information. The spatially contextual MRF model is time- 
consuming and cannot describe long-distance interactions between pixels, i.e., a macro texture model. The 
region-based MRF, on the other hand, can use region-based information to capture the macro texture pattern, 
but it suffers from irregular spatial context. To achieve accurate measurement of the interface in real-time, we 
use an approach based on Gaussian mixture models and MRF based on unsupervised image segmentation 
[15]-[17]. Image processing tasks such as image segmentation have proven to be suitable for MRFs, [18]. 
One for the background and one for the foreground GMM, replaced the user's seed histograms and to fix the 
transparency on the edges of segmented objects, a border matting algorithm was developed [19]-[21]. 

The region-level MRF model has been proposed with an approach for the fusion using the FCM 
clustering algorithm based on the MRF, multispectral and panchromatic images are exploited for local spatial 
information [22]-[25]. Their concept of label uncertainty, however, cannot be applied to generic MRF 
models. A general concept of uncertainty in the context of image segmentation can be defined as the variance 
or inequality [26] in the label assignments resulting from the label image's distribution posterior. To solve the 
MAP problem, the iterative conditional model (ICM) is used; this method is compared to traditional 
expectation-maximization (EM) and MRF image segmentation techniques using synthetic and real images 
[27], [28]. The approximation algorithms they induce are iterative and very computationally intensive. The 
joint field of expectation of maximization field and MRF are designed to optimize image segmentation 
results, which can improve image segmentation accuracy. 

The rest of the research paper is organized as shown in: Section 2 provides a brief overview of the 
MRF model, while section 3 presents the new regional feature of the MRF model as well as details of the 
EMRF model. Section 5 describes the algorithm. Section 3 contains the results of the experiments, and 
section 4 contains the conclusion. 


2. MATERIALS AND METHOD 
2.1. Framework MRF 

Markov random field theory provides a handy and coherent framework method for modeling image 
pixels and linked features are examples of context-dependent entities. Improved maximum a posteriori 
expectation-maximization (MAP-EM) image segmentation algorithm that calculates the proportion of tissue 
mixture within each picture voxel as well as tissue distribution statistical model parameters [29]. It may be 
viewed as a labeling issue for segmentation, with the segmented result generated by maximizing the posterior 
estimation of the labeling problem. By synthesizing the maximum, a posteriori (MAP) solution and 
incorporating prior information, the Bayesian framework aids in the generation of statistical inferences. The 
purpose of the MAP solution is to maximize the posterior probability, which is written as: 

MRF framework: Given a group of pixels S = {s, ...Sm} and group of labels A= {], ...1,} and 
neighborhood system N, find the mapping of S to A. 

Let K be the configuration for labels, 


M = {1, ...1,}, mjeA is the label for s;. 


M is MRF concerning N iff: 
— Positivity: PCM = m) > 0VmeM. 
—  Markovianity: 


_m fs = _m os 
P(M, = /M, =m,,Vr# s) = P(M, = /M, = m,, Vre N) (1) 
Let G be the observed image. 
- G=q(H(M) +N) 


Where H=camera transfer function and @ =Recorder distortion, h is assumed to be LSI and @ is invertible 


nonlinear function and N is additive noise. In the framework of restoration: Given G what is F?; P (m = 7 a 


2), Maximum likelihood of: M = m given; G = g. From Bays rule: 


Markov random field model and expectation of maximization for images ... (Lalaoui Lahouaoui) 


774 0 ISSN: 2502-4752 


P(M===g)aP(G=£=M)P(M=m) (2) 


where P (G=£ =M)= data model, P(M = m)= prior and P (m= = g) = Aposteriori distribution. Need 


to maximize posterior (MAP) distribution: 

P(k) = (1/2Z) exp (— UC m)/T) (3) 
UCM) = YicenV-(m) IS energy function, where V,: clique Potentiel: 

Z = Yir(1/Z) exp (— U(M)/T) (4) 


Z: partition Function, and T= Temperature. 
Hammersley Clifford theorem M is MRF on S concerning N if and only if M is Gibbs random 
field on S forl. Relates conditional distribution (Local characteristic joint distribution (Gibbs measure)). 


argmaxm € M P(f|d) (5) 


If computed over all cliques, it would take a long time. In practice, it is common to execute the computation 
over single and pair site cliques. 


2.2. The MRF model 

The Markov field technique attempts to account for term dependency. It's also known as an 
undirected graphical model. This method is widely used in machine learning. This method is used in IR to 
describe term dependencies, such as dependence sequence and complete dependency. 


2.3. The proposed algorithms 

In the EM algorithm simulated field. Input: an observation y, some iterations T. Output: an estimate 
of the likelihood maximum estimator @MLE. Initialization: start to form an initial guess@ = (8, gp), do 
Neighbourhood restoration: draw 6 from n(.| y, 0%", @). 

E-step: Estimate the expected value for each latent variable. We compute: 


B,D) = LicoDxy P(x RN Yor WO, A loggy(yilxi) (6) 
G1 (O) = Lier Ux; P(x RN Yu WO, L) logP (xi; RQ2y,Yi ®, L) (7) 


M-Step. Optimize the parameters of the distribution using maximum likelihood. set of parameters, k denotes 
the Gaussian component of the mixture model, K denotes the total number of Gaussian distributions in the 
mixture model. Wm is the weight of each Gaussian distribution that meets M, m=1, and wn=1. 

6: old, 0’ : New, 


Porty,x) 
96,8") = Dy Per) rogpoers) = LY Gory BPO) 


=D 0’ (2) log (po®)Pew) = Ly Po'®) [logpe (2) + logpe(s)| 
= YyPy@) logps (2) + Ly py® logpo(x) (8) 


logpe (¢ \st eeu 4: logpe(x) <x-1 


06,8") < Dy per (2) (2-1) + Lype (Z)@-D (9) 
08,6") < Dy pe (2).2— Ly per (2) + Ly per (%)-x— Ly per (2) (10) 
06,6") < Dypor (2).2+ Dy por (%).x- 2 Dy per (2) (11) 
0(8 6") < Dy por (%).2 + (— 2) Ly por (2) (12) 


end, 
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a = (0M, gM) (13) 
that is, 
P(X = x) = (1/Z)exp(—U(QX)) (14) 
Where Z = YxU(x) (5) 


is the normalization factor and, 

U(x) = Ys € SU(xs, xNs) (16) 
is the energy function. Here: 

U(xs, xNs) = Vien V(xs, xt) (17) 
where Ns is a group of sites adjacent to site s and V(xs,xt) is the potential function between site s and site t, t 
Ns. Typically, the multilevel logistic (MLL) model [30] is used to define the potential function: V(xs,xt) 
which is: 


V(xs,xt) = Bifxs # xt = —BPifxs + xt, (18) 


where B >0 denotes the potential parameter and t€ Ns. According to the MLL model, P(X) has a large value 
if the local neighbor labels are the same, otherwise it has a small value It is used to calculate the probability 
that the observed image Y corresponds to the given realization X = x. That is: 


P(Y[X = x) = [I sesP(¥s|X) = [I sesP(Yslxs) (19) 


segmentation consists in estimating 1. knowing the observed image a. You can use a Bayesian approach 
which consists in finding the most probable class in the sense of the probability at posterior obtained by 
maximizing P(ys = p|xs =I). It is therefore a matter of constructing P(xs = I'| ys = p), then fix a prior law 
on p, denoted P(u =I), from which we calculate P(ys = p|xs = I) by of Bayes. 


P(ys|xs) = (2m) — D/2T — 1/2exp — 1 2(ys — pw) TT — 1(ys — p) (20) 


3. RESULTS AND DISCUSSION 

The obtained results are presented in the Figure 1 to Figure 3, where we have presented the 
segmentation evaluation criteria in Table 1. It is concluded that the ICM algorithm based on the EM method 
performs better compared to the ICM based on the method for K-means ima+ge real. In this section, the 
proposed algorithm is applied to synthetic and real images. We present the results of the application of 
improved MRF based on modified EM algorithms. Our proposed technique has been tested on several 
different images. Figure 1 shows the input images and image segmentation obtained by our method. It is 
clear that the results using the proposed technique had good quality and were better than those using the 
classical method. The performance of these algorithms is compared to the standard versions of the MRF and 
EM algorithms. We present the results obtained for image segmentation by Markov fields based on EM and 
changed EM algorithms. For output image evaluation, not only visual comparison but also quantitative 
measurements such as entropy, mean square error (MSE), and peak signal to noise ratio (PSNR) are used to 
evaluate the output image. To make an objective comparison of the different proposed methods, we propose 
to use the following as evaluation criteria: mean squared error (MSE), signal to noise ratio PSNR, normalized 
cross-correlation (NCC), structural content (SC) and normalized absolute error (NAE). We present in tables 
the values of the validation criteria for all the proposed methods. This will be followed by an interpretation of 
the results to be able to make a comparison between the used segmentation methods. First, we define the 
evaluation criteria that we mentioned above. The methods of assessing image quality can be classified into 
three broad categories: The most widely used quantitative measures are the mean squared error (MSE), the 
peak signal-to-noise ratio (PSNR), and the signal-to-noise ratio (hat the ICM algorithm based on the EM 
method performs better compared to the ICM algorithm based on the method for K-means image real. In this 
section, the proposed algorithm is applied to synthetic and real images. We present the results of the 
application of improved MRF based on modified EM algorithms. Our proposed technique has been tested on 
several different images. Figure 1 shows the input images and image segmentation obtained by our method. It 
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is clear that the results using the proposed technique had good quality and were better than those using the 
classical method. The performance of these algorithms is compared to the standard versions of the MRF and 
EM algorithms. We present the results obtained for image segmentation by Markov fields based on EM and 
changed EM algorithms. For output image evaluation, not only visual comparison but also quantitative 
measurements such as entropy, mean square error (MSE), and PSNR are used to evaluate the output image. 
To make an objective comparison of the different proposed methods, we propose to use the following as 
evaluation criteria: MSE, PSNR, NCC, SC, and NAE. It is based on measuring the MSE computed between 
the original pixels and the degraded ones: The MSE is equal to 0 if a model has no errors. Its value rises in 
proportion to the model error. Additionally known as the mean squared deviation (MSD), the mean squared 
error. 

Where (M x N) is the size of the image, Ip and Ip are respectively the amplitudes of the pixels on the 
original and degraded images. The eye likely takes more account of errors in large amplitudes, which 
promotes the quadratic measure. The PSNR is to quantify the performance of the algorithms by measuring 
the quality of reinsertion of the segmented image compared to the original image. 

PSNR levels for high-quality images typically range between 30 and 40 dB. A greater SC rating 
indicates that the image is of low quality. The NCC calculation compares the processed picture to the 
reference image. Finally, the greater the value of the maximum difference, the better the image quality. This 
quality metric may be stated as follows: This measure of quality can be stated as follows: A higher price 
NAB demonstrates that the image quality is low. 


3.1. Synthetic images segmentation 

In this part, we started to apply the two segmentation methods to the synthetic images, and after that; 
we evaluated these methods by evaluation criteria that we saw previously. From the values in Table 1, 
(PSNR=7.9280) and (SC=1.5034) we see that the quality of the segmented image is poor by the ICM method 
based on the EM algorithm versus ICM. An infinite PSNR value corresponds to an undegraded image. In 
addition, this value decreases according to the degradation of the PSNR and therefore links with MSE the 
maximum energy of the image. Regardless of the value of MSE, the value of PNSR is large. In Figure 1 (a), 
we presented the synthetic image original, these results obtained after segmentation with MRF-based EM (in 
Figure 1(b), and Segmentation with MRF-based EM algorithm in Figure 1(c). From the values obtuned in 
Table 1 (PSNR=7.9280) and (SC=0.2115) where we can see that the quality of the segmented image is poor 
from the ICM method based on the EM algorithm. Table 1 presented the two algorithms used for synthetic 
image segmentation with EM and EM modified. 


Figure |. Result obtuned of simulation methods for synthetic images, (a) image original, (b) segmentation 
with MRF-based EM, and (c) segmentation par MRF-based modified EM algorithm 


Whatever the value of MSE, the value of PNSR is great. In Table 1, we present also the values of 
different criteria used here such as MSE, NCC, SC, PSNR, and NAE for the EM and EM modified used in 
Markov field random. Figure 2 presented the bare for criteria MSE, PSNR, NCC, SC and NAE obtained from 
the synthetic image segmentation algorithm. 


Table 1. Values of the evaluation criteria for image segmentation 
Criteria _ EM modified EM 


MSE 1.6607 1.1209 
PSNR 7.9280 7.6352 
NCC 2.0332 1.9389 

SC 0.2115 0.2501 
NAE 1.0332 0.9389 
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modified EM and EM 
GEES Modified EM 
EM 


The values of criteria 


2} I | 
od | —— Mn | 
MSE PSNR NCC sc NAE 
Different criteria such as MSE, NCC, SC, PSNR, and NAE 


Figure 2. We present here the accuracy criteria such as (MSE, PSNR, NCC, SC, and NAE) obtained with 
MRF based on EM, and MRF based on the modified EM for the image synthetic 


3.2. Real image segmentation 

In this part, we started to apply the two segmentation methods to real images and synthetic images 
and then we evaluated these methods by evaluation criteria that we have seen previously. The obtained 
results are depicted in Table 2 and it shows that the ICM algorithm based on the modified EM method is 
better compared to ICM based method for EM on the actual images. To benefit from the advantages of the 
two methods presented in the previous paper, we decided to use them to segment the images, and we made a 
table for each result. Table 2 presents the values of some evaluation criteria on the segmentation results with 
the two optimal values (taking into account the direction of variation of the criteria). 

Table 2 presented the value of different criteria used here such as MSE, NCC, SC, PSNR, and NAE 
for the EM and EM modified used in Markov field random. Figure 3 presented the two algorithms used for 
synthetic image segmentation with EM and EM modified. In Figure 3(a), we presented the reel image 
original, these results obtained after segmentation with MRF-based EM in Figure 3(b), and Segmentation 
with MRF-based modified EM algorithm in Figure 3(c). From the values in the Table 2 (PSNR = 6.925) and 
(SC = 0.2349) in can see that the quality of the segmented image is poor by the ICM method based on the 
EM algorithm versus ICM. 

Figure 4 presents the bare for MSE, PSNR, NCC, SC, and NAE criteria obtained from the synthetic 
image segmentation algorithm. To benefit from the advantages of the two methods presented in the previous 
section, we decided to use them to segment the images. And we made a table for each result; this table shows 
the values of some evaluation criteria on segmentation results both optimal values (considering the direction 
of change of criteria). We can see in Figure 4 that the interest of the segmentation method KMeans it will 
keep precise borders. The outline image as it appears is a typical image for this kind of segmentation with 
similar EMs. The values of KMeans are large relative to the value of EM. Therefore, the quality is good and 
the best EM algorithm. It is possible to deduce that the ICM algorithm based on EM method is better 
compared to ICM based on the method for K-means images the values of the criteria give the information 
that the EM algorithm is better than the K-means algorithm. 


Table 2. The Values of criteria for image segmentation 


Criteria EM modified EM 
MSE 1.3199 1.5917 
PSNR 6.925 6.1122 
NCC 1.9689 1.8724 

SC 0.2349 0.2305 
NAE 0.96890 0.8724 


(a) (b) 


Figure 3. Result obtained with simulation for real images, (a) image original, (b) segmentation with MRF 
based on EM, and (c) MRF based on modified EM 
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Figure 4. We present here the accuracy criteria such as (MSE, PSNR, NCC, SC, and NAE) obtained with 
MRF based on EM, and MRF based on the modified EM for the image real 


4. CONCLUSION 

In current work, we have considered various segmentation tasks in an image labeling framework, 
where the challenge is simplified to labeling pixels. The labels were modeled using MRF in a probabilistic 
approach, and inference is performed using Bayesian estimation, specifically a MAP. MRF modeling has the 
advantage of allowing a priori information to be "coded" locally by click potentials. We exploit the results 
obtained by different methods such as Markov fields based on the EM and modified EM algorithms for 
image segmentation. Experimental evaluations on multiple reference sources and target datasets indicate that 
this is efficient for image segmentation. We proposed a method for improving the image's grayscale and 
image detail by using quantitative measures to analyze the output image, such as entropy, MSE, and PSNR. 
Image processing is a research subject located between computer science and signal processing. In current 
work, we have presented various algorithms for based-region segmentation. Whether sequential or mutual, 
cooperative segmentation uses a single type of information per region to allow better consideration of the 
characteristics of the objects in the image. The region-wise approach consists in bringing together the 
associated pixels in a homogeneous region; it is based on the notion of similarity. It's quick and easy, but 
using only local MRF information based on EM to assess these results, we used the following criteria: MSE, 
PNSR, NCC, SC, MD, and NAE, and two algorithms: EM and EM amended. In this work, we presented the 
results of the region segmentation algorithm based on Markov fields. The obtained results show that the MRF 
method based on modified EM is better than that based on classical EM. Nevertheless, many image 
segmentation tasks can be formulated under these constraints. 
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